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Figure  Captions 


Figure  1.  Flowchart  showing  the  major  elements  of  the  ionospheric  model. 

Figure  2.  PIIE  observations  of  electron  density  fluctuations,  total  electron 
density,  electron  temperature,  and  energy  flux  of  precipitating  electrons  vs. 
time  (seconds  after  launch),  altitude,  and  UT.  Vertical  line  indicates  time  of 
the  model  calculation  Inside  the  arc  (284  SAL). 

Figure  3.  Differential  flux  spectrum  of  downgoing  electrons  input  to  the  model 
corresponding  to  the  time  indicated  by  the  vertical  line  in  Figure  2. 

Figure  4.  Example  of  typical  electron  flux  observed  outside  arc. 

Figure  5.  Model  calculations  of  total  electron  density,  CH-  density,  and 
electron  temperature  resulting  from  the  flux  in  Figure  3. 

Figure  6.  Hall  and  Pedersen  electric  conductivities  resulting  from  the  flux  in 
Figure  3. 

Figure  7.  Relative  densities  of  N+,  N0+,  and  0+  observed  during  the  flight. 

Figure  8.  Relative  densities  of  N2+  and  O2+  observed  during  the  flight. 

Figure  9.  Calculated  electron  density  profiles  corresponding  to  increasing  flux 
duration  (c.f.  Figure  3). 

Figure  10.  Same  as  Figure  9  except  starting  from  a  lower  initial  electron 
density. 

Figure  11.  Decay  of  electron  density  profile  after  flux  in  Figure  2  is  turned 
off. 

Figure  12.  Development  of  electron  density  profile  using  electron  flux  outside 
arc  (c.f.  Figure  4). 

Figure  13.  Same  as  Figure  12  except  starting  from  a  lower  initial  electron 
density. 

Figure  14.  Calculated  electron  density  profiles  using  flux  inside  arc  and  a 
cold  neutral  atmosphere  (Tjj^ySOK). 

Figure  15.  Same  as  Figure  13  except  with  a  hotter  neutral  atmosphere  (Tjj“1200K). 
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1.  Introduction 

A  rocket  experiment  conducted  in  Greenland  on  March  15,  1985,  the  Polar 
Ionospheric  Irregularities  Experiment  (PIIE),  took  detailed  measurements  of 
important  quantities  as  it  passed  over  an  F-layer  polar  cap  auroral  arc. 
Descriptions  of  the  instruments  and  observations  during  the  flight  can  be 
found  in  Weber  et  al.  (1987).  Data  from  the  mission  provided  key  inputs  and 
valuable  checks  for  an  ionospheric  model  capable  of  handling  soft  electron 
fluxes  often  seen  over  polar  cap  auroral  arcs.  The  original  model  codes  for 
differential  electron  flux  throughout  an  ionospheric  column  (flux  transport) 
were  developed  by  Strickland  et  al.  (1976);  local  chemistry  and  t>+  diffusion 
codes  were  added  to  find  ion  densities  and  optical  emissions.  This  model  was 
used  by  Weber  et  al.  (1985)  to  identify  the  source  of  F-region  plasma 
enhancements  at  tlie  equatorward  boundary  of  the  auroral  zone.  Later,  we 
included  an  electron  temperature  calculation  with  heat  conduction  along 
magnetic  field  lines,  permitting  us  to  more  accurately  decermine  temperature 
dependent  parameters  such  as  chemical  rate  and  OP  diffusion  coefficents. 

The  flowchart  in  Figure  1  shows  the  principal  features  of  the  model.  In 
its  present  setup  the  code  has  two  distinct  parts:  part  1  in  the  upper  section 
(transport)  calculates  production  rates  of  important  ion  and  neutrals  that 
result  from  a  flux  of  incident  electrons  at  an  upper  altitude  boundary;  part  2 
uses  the  production  rates  to  calculate  ionospheric  densities,  temperature,  and 
electric  conductivities.  The  flux  transport  program  requires  a  neutral 
atmosphere  model  consisting  of  N2,  O2,  and  0  and  a  table  of  cross  sections 
describing  important  inelastic  scattering  processes  involving  electrons  and 
the  ambient  neutrals. 

Strickland  et  al.  (1976)  demonstrate  how  the  Boltzmann  equation  is 
reduced  to  a  one-dimensional,  steady  state,  transport  equation  for 


differential  electron  flux  over  altitude.  Scattering  terms  involve  both 
elastic  and  inelastic  collisions.  Since  electrons  cannot  gain  energy  in  the 
model,  no  external  forces  are  considered;  secondaries  are  produced 
isotropically  and  contribute  to  the  thermal  population.  Because  the  flux  at  a 
given  energy  and  pitch  angle  can  result  only  from  electrons  at  tlie  same  or 
higher  energies,  the  code  steps  successively  to  lower  energies  in  one  pass. 

A  matrix  representation  of  the  transport  equation  is  solved  for  2U  pitch 
angles  over  an  energy  grid  appropriate  for  the  incident  differential  flux 
(el/cm-s-sr-ev) .  An  exact  eigenvalue  solution  is  found  at  some  intermediate 
altitude,  which  leads  to  solutions  at  other  altitudes  from  standard 
perturbation  techniques.  The  code  is  relatively  fast  yet  more  sophisticated 
than  empirical  range  theoretic  models  (Rees ,  1963).  It  does  not  require  a 
matching  of  solutions  between  high  and  low  energy  ranges  as  do  schemes 
employing  a  Kokker-Flanck  approximation  at  high  energies  (Banks  et  al.,  1974). 
Mantas  and  Walker  (1975)  use  a  similar  approach  that  differs  somewhat  in  the 
representation  of  flux  vs.  energy.  An  energy  conservation  test  compares  the 
energy  lost  by  the  electron  flux  against  the  energy  expended  in  ionization  and 
excitation  of  the  neutral  constituents. 

Chemical  production  rates  and  the  electron  plasma  heating  rate  calculated 
by  the  flux  transport  code  become  inputs  to  the  ionospheric  chemistry  and 
temperature  codes  in  the  bottom  of  Figure  1.  These  inputs  are  constant  in 
time  but  can  be  switched  off  to  simulate  decay  in  the  absence  of  an  ionizing 
flux.  Starting  from  initial  values  for  the  electron  temperature  and  densities 
of  various  neutral  ionic  species,  the  main  chemistry  code  integrates  the 
chemical  rate  equations  using  a  predictor-corrector  method  with  an  adaptive 
time  step.  Because  0+  is  the  only  ion  appreciably  affected  by  diffusion  along 
magnetic  field  lines  it  is  handled  separately  in  concert  with  the  other 


species.  The  calculated  electron  density  (total  ion  density)  then  determines 
the  electron  heat  conduction  coefficient  and  the  total  electron  plasma  heating 
rate . 

A  Crank-Nicolson  finite  differencing  scheme  solves  both  the  heat 
conduction  and  diffusion  equations;  the  method  is  implicit,  meaning  that  a 
solution  involves  unknown  values  at  adjacent  altitudes.  Solutions  are  always 
stable  (Gerald  and  Wheatley,  1984)  but  at  the  price  of  having  to  solve  a  si-.t 
of  simultaneous  equations.  The  diffusion  equation  is  relatively 
straightforward  since  the  coefficients  depend  only  on  temperature  and  other 
parameters  which  are  updated  at  regular  time  intervals.  The  temperature, 
however,  requires  the  solution  of  a  non-linear  equation  with  coefficients  that 
depend  on  the  temperature.  With  short  time  steps  the  problem  can  be  managed 
by  assuring  that  the  coefficients  change  slowly  enough  to  be  linearly 
projected  to  the  end  of  each  time  step.  Alternatively,  other  researchers  have 
used  Newton's  method  to  solve  a  different  formulation  of  the  non-linear, 
diffusion  type  partial  differential  equation  (Hastings  and  Roble,  1977; 

Mantas  et  al.,  1981). 

The  resulting  Tg  is  then  used  to  recalculate  temperature  dependent  rate 
coefficients  and  the  0+  diffusion  coefficient.  The  ion  temperature,  Tj^,  is 
not  calculated  now  by  the  code  but  is  set  between  the  electron  and  neutral 
temperatures.  Our  experience  indicates  that  a  loop  through  the  chemistry, 
diffusion,  and  temperature  codes  every  second  for  the  first  2  minutes  in 
ionosphere  time  yields  reasonable  results-  After  2  minutes,  when  quantities 
change  more  slowly,  the  update  time  is  lengthened  to  10  seconds.  Pederson  and 
Hall  electric  conductivities  are  also  calculated  from  expressions  found  in 
Val lence-Jones  (1974)  with  the  derived  electron  temperature  and  density 


The  present  model  is  limited  to  one  dimension  along  a  magnetic  field  line 
and  neglects  horizontal  transport.  This  is  adequate  to  describe  an  auroral 
arc  several  kilometers  wide  that  is  either  stationary  or  has  its  production 
source  moving  with  the  convection  velocity.  Real  physical  conditions  are  more 
complicated  than  these  simple  assumptions;  the  prime  aim  of  this  work  is  to 
find  out  how  well  the  simplified  model  fares  in  a  carefully  documented  case 
study.  Other  ionospheric  models  with  rather  complex  chemistry  have  been 
developed  over  the  past  15  years  (e.g.,  Jones  and  Rees ,  1963;  Roble  and  Rees, 
1977;  Schunk  et  al.,  1986).  The  main  advantage  of  this  model  centers  on  the 
detailed  specification  of  production  rates  resulting  from  a  measured 
differential  electron  flux.  Global  or  limited  two  dimensional  models  would 
require  many  such  calculations  in  addition  to  Horizontal  transport  in  the 
chemistry  model.  A  parameterization  of  the  flux  transport  code  to  handle  many 
situations  is  the  next  logical  step  in  the  development  of  this  particular 
model . 

2.  PIIE  Measurements  and  Model  Results 

Figure  2  (Figure  5  from  Weber  et  al.  (1987))  displays  measurements 
taken  by  the  instrumented  rocket  payload  as  it  passed  over  a  complicated 
system  of  F-layer  auroral  arcs  near  apogee.  The  vertical  line  indicates  the 
time  chosen  for  model  comparisons.  The  electron  density  had  a  local  maximum 
of  3  X  10^  el/cm^  and  Tg  was  about  3000K  at  415  km  altitude.  Electrons 
precipitating  into  the  arc  carried  an  energy  flux  of  about  .1  erg/cm^s,  nearly 
an  order  of  magnitude  greater  than  in  adjacent  regions.  The  differential 
electron  flux  spectrum  at  this  time  appears  in  Figure  3.  Although  only  a  .5  s 
average  sample,  the  spectrum  is  representative  of  a  5  s  interval  near  tne 
center  of  the  arc.  During  the  time  wJien  the  rocket  was  travelling  almost 
horizontally  near  apogee,  one  of  the  electrostatic  analyzers  sampled  pitch 
angles  near  45°.  Since  there  were  no  measurements  taken  simultaneously  at 


other  pitch  angles  less  than  90°  we  assumed  an  isotropic  distribution  of 
downgoing  electrons.  Outside  the  arc  (Figure  4)  fluxes  were  2  orders  of 
magnitude  lower  but  decreased  more  gradually  at  higher  energies. 

Model  results  inside  the  arc  (Figure  5)  indicate  a  peak  electron  density 
of  about  5  X  10^  el/cm^  at  300  km  altitude  and  3.5  x  10^  el/cm^  at  the  rocket 
apogee  altitude;  Tg  at  145  Xm  altitude  was  calculated  to  be  slightly  greater 
than  20U0K.  0+  ions  dominate  at  high  altitudes  but  are  replaced  by  increasing 
amounts  of  NtH-,  O2+,  N+,  and  N2+  at  lower  altitudes.  These  results  start  froia 
an  Initial  n^  =  10^  el/cm^  at  all  altitudes  and  Tg  =  3000K.  at  600  km  altitude, 
after  which  production  by  the  precipitating  electrons  is  turned  on  for  2 
minutes.  More  realistic  initial  densities  estimated  from  the  upleg  density 
profile  outside  the  arc  yield  somewhat  higher  densities  (~5xl0^  el/cm^).  The 
calculated  Tg ,  however,  is  not  greatly  affected  when  appreciably  different 
initial  conditions  are  used.  Hall  and  Pedersen  electric  conductivities 
calculated  inside  the  arc  are  shown  in  Figure  6.  Even  though  most  ionization 
occurs  above  200  km  both  conductivities  peak  low  in  the  E  region  where 
collisions  with  neutrals  rapidly  increase. 

Although  the  model  electron  density  compares  reasonably  well  with  the 
data,  the  electron  temperature  is  somewhat  lower  than  measured.  Because  the 
only  source  now  incorporated  in  the  model  is  collisional  heating  by  energetic 
electrons,  this  discrepency  may  in  part  stem  from  other  types  of  heating. 

More  fundamentally,  heat  conduction  is  only  part  of  the  full  energy  balance 
equation  (Schunk  and  Walker;  1970>;  terms  representing  other  forms  of  heat 
transport  could  modify  the  results.  For  example,  Schunk  et  al.  (1987) 
demonstrate  that  thermoelectric  cooling  by  field-aligned  currents  becomes 
significant  only  for  rather  large  return  currents  of  thermal  electrons  on  the 
order  of  lU  pA/tn2.  Weber  et  al.  (19871  estimate  maximum  field-aligned 


currents  of  -'3  yA/mi!  within  the  arc,  suggesting  a  minimal  contribution  from 
this  effect. 

Relative  densities  of  several  ionic  species  can  also  be  compared  with 
model  results.  Figures  7  and  8  show  ion  currents  for  several  species  measured 
by  the  ion  mass  spectrometer.  Since  the  current  is  proportional  to  the  ion 
density,  relative  densities  of  the  species  can  be  computed.  Density  ratios  of 
0+/N+  and  N+/NO+  were  about  10  to  1  and  3  to  1,  respectively.  With  an  initial 
profile  containing  only  only  0+  (=  1x10^  cm“^  at  all  altitudes)  the 
corresponding  model  ratios  after  2  min  are  close  to  the  measured  ones.  With  a 
more  realistic  initial  profile  (0+  =  1.13  x  10^,  N0+  =  1.42  >;  Id^,  U2+  = 

7.41  X  lU'^,  N+  =  3.43  X  lo^,  and  N2+  =  9.59  x  10°  cm~^,  all  at  AOO  km), 

0+/h+  =  12  and  N+/N0+  =1.2  after  2  min.  While  the  total  ion  density  and 
0+/N+  ratio  at  4U0  km  are  relatively  insensitive  to  initial  conditions,  otner 
minor  species  show  greater  dependence  on  initial  densities,  especially  several 
minutes  after  t  =  0.  For  exaiaple ,  Figure  8  indicates  that  O2+  dominates  N2+ 
at  400  km  while  tiie  model  has  N2+  =  1.02  x  10^  and  O2+  =  5.48  x  lu^  aa“^  when 
0+  is  the  only  species  initially  present.  The  realistic  initial  condition 
leads  to  N2+  =  9.51  x  10^  and  U2+  =  6.26  x  10^  cm“^,  a  much  closer 
correspondence . 

3.  Discussion 

Densities  associated  with  an  arc  lasting  a  few  minutes  to  several  tens  of 
minutes  depend  both  on  the  initial  ionosphere  and  the  duration  of  enhanced 
production.  Although  direct  procuction  by  precipitating  electrons  leads  to  a 
rapid  buildup  ot  ionization,  equilibrium  resulting  from  chemistry  and 
diffusion  is  reached  slowly  (several  tens  of  minutes  in  the  K-region  vs. 


minutes  in  the  E-region).  Starting  with  rather  artificial  initial  conditions, 
the  next  two  examples  (Figures  9  and  10)  illustrate  how  the  electron  density 
profile  develops  when  a  flux  like  that  observed  during  FtlE  is  turned  on 
continuously  for  30  minutes. 

In  Figure  9,  0+  ions  alone  are  present  initially  with  a  constant  height 
profile  of  10^  cm“^.  From  120  to  200  km  altitude  production  and  loss  quickly 
balance  and  equilibrium  densities  are  reached  after  10  minutes  or  so.  Below 
120  km  production  is  slight  and  the  density  actually  decreases.  Ion 
production  dominates  above  200  km,  causing  a  steady  increase  in  electron 
density  that  is  limited  by  CH-  diffusion  until  equilibrium  is  reached  several 
hours  later.  Starting  with  an  0+  density  of  10^  cm“^.  Figure  10  indicates 
that  the  final  density  profile  after  30  minutes  is  not  different  than  seen  in 
Figure  8,  except  at  the  lowest  altitudes.  From  130  to  200  km  the  equilibrium 
densities  are  nearly  equal;  above  200  km  the  profiles  are  similar  but  lag 
behind  those  in  Figure  8  because  of  the  smaller  initial  density.  However,  on 
time  scales  longer  than  5  minutes,  the  resulting  densities  above  100  km  agree 
to  within  20%  even  when  drastically  different  initial  conditions  are  assumed. 

Figure  11  shows  what  happens  it  the  flux  is  turned  on  for  15  minutes  and 
then  suddenly  switched  off,  after  which  direct  production  ceases  and 
ionization  decays,  subject  only  to  local  chemistry  and  diffusion.  Following 
an  abrupt  drop  caused  by  rapid  losses,  the  electron  density  slowly  decays, 
especially  at  and  above  the  altitude  peak.  Here  diffusion  flattens  the 
profile  and  maintains  the  density  above  400  km  where  losses  by  recombination 
are  minimal. 

For  comparison.  Figures  12  and  13  show  electron  density  profiles 
generated  by  fluxes  measured  outside  the  main  auroral  arc  at  about  306  sec 
into  the  flight  (Figure  4).  Initial  0+  densities  are  10^  and  10^  cm“^  as 


before.  Here  the  initial  density  conditions  clearly  affect  the  profiles  even 
after  several  minutes  have  elapsed.  Only  in  the  E-region  where  local 
chemistry  is  relatively  fast  do  the  densities  approach  the  same  equilibrium, 
and  neither  case  approximates  the  densities  observed  in  the  downleg  segment. 
Since  the  electron  energy  flux  measured  from  420  k.m  to  290  km  nowhere 
approaches  the  more  energetic  flux  seen  earlier,  it  seems  likely  that  the 
observed  electron  density  results  from  a  combination  of  ionization  remaining 
after  local  production  ceased  and  that  brought  in  by  convection. 

A  future  goal  of  this  study  is  to  find  the  relationship  between  large 
scale  features  in  the  polar  cap  ionosphere  and  scintillations  of  satellite  to 
radio  transmissions.  Scintillations  are  caused  by  plasma  density 
irregularities  with  scale  sizes  from  100  m  to  10  km  (Weber  and  bachau,  19bb). 
These  irregularities  have  been  associated  with  plasma  instabilities  generated 
by  density  gradients  in  enhanced  plasma  regions  (blobs)  hundreds  of  km  across 
or  with  velocity  shears  (Basu  et  al.,  1966).  /lence,  density  features 
resulting  from  direct  production  or  transport  play  an  important  role  in  the 
problem.  For  example,  observed  changes  in  scintillation  behavior  with  the 
solar  cycle  (Aarons ,  1982)  may  be  connected  with  the  way  these  features  are 
generated.  A  model  can  probe  how  similar  fluxes  affect  neutral  atmospheres 
with  different  exospheric  temperatures.  In  Figure  14,  a  cold  atmosphere  like 
that  found  during  solar  minimum  produces  a  slightly  higher  density  at  lower 
altitudes  (cf.  Figure  9,  =  lOOOK).  A  hot  neutral  atmosphere  has  peak 

dens' ty  above  3U0  km  and  a  more  gradual  decrease  toward  higher  altitudes 
(Figure  15).  Differences  between  profiles  may  help  identify  which  processes 
generate  irreguiatities  and  the  scale  sizes  likely  to  be  observed. 
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